Molecular dynamics simulations of D2O ice photodesorption 

C. Arasa 

Leiden Observatory, Leiden University, P. 0. Box 9513, 
2300 RA Leiden, The Netherlands and 
Gorlaeus Laboratories, Leiden Institute of Chemistry, 
Leiden University, P. 0. Box 9502, 2300 RA Leiden, The Netherlands 

S. Andersson 

SINTEF Materials and Chemistry, P.O. Box 4760, 7465 Trondheim, Norway and 
Department of Chemistry, Physical Chemistry, 
University of Cothenburg, 41296 Gothenburg, Sweden 

H. M. Cuppen and E. F. van Dishoeck 

Leiden Observatory, Leiden University, 
P. 0. Box 9513, 2300 RA Leiden, The Netherlands 

G. J. Kroes 

Gorlaeus Laboratories, Leiden Institute of Chemistry, 
Leiden University, P. 0. Box 9502, 2300 RA Leiden, The Netherlands 

(Dated: January 20, 2013) 



1 



Abstract 

Molecular dynamics (MD) calculations have been performed to study the ultraviolet (UV) pho- 
todissociation of D2O in an amorphous D2O ice surface at 10, 20, 60, and 90 K, in order to 
investigate the influence of isotope effects on the photodesorption processes. As for H2O, the 
main processes after UV photodissociation are trapping and desorption of either fragments or D2O 
molecules. Trapping mainly takes place in the deeper monolayers of the ice, whereas desorption oc- 
curs in the uppermost layers. There are three desorption processes: D atom, OD radical, and D2O 
molecule photodesorption. D2O desorption takes places either by direct desorption of a recombined 
D2O molecule, or when an energetic D atom produced by photodissociation kicks a surrounding 
D2O molecule out of the surface by transfering part of its momentum. Desorption probabilities 
are calculated for photoexcitation of D2O in the top four monolayers and compared quantitatively 
with those for H2O obtained from previous MD simulations of UV photodissociation of amorphous 
water ice at different ice temperatures [Arasa et al., J. Chem. Phys. 132, 184510 (2010)]. The 
main conclusions are the same, but the average D atom photodesorption probability is smaller than 
that of the H atom (by about a factor of 0.9) because D has lower kinetic energy than H, whereas 
the average OD radical photodesorption probability is larger than that of OH (by about a factor 
of 2.5-2.9 depending on ice temperature) because OD has higher translational energy than OH for 
every ice temperature studied. The average D2O photodesorption probability is larger than that of 
H2O (by about a factor of 1.4-2.3 depending on ice temperature), and this is entirely due to a larger 
contribution of the D2O kick-out mechanism. This is an isotope effect: the kick-out mechanism is 
more efficient for D2O ice, because the D atom formed after D2O photodissociation has a larger 
momentum than photogenerated H atoms from H2O, and D transfers momentum more easily to 
D2O than H to H2O. The total (OD -|- D2O) yield has been compared with experiments and the 
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total (OH + H2O) yield from previous simulations. We find better agreement when we compare 
experimental yields with calculated yields for D2O ice than when we compare with calculated yields 
for H2O ice. 

PACS numbers: 
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I. INTRODUCTION 



The formation of molecules in the interstellar medium (ISM) can proceed through several 
kinds of reactions. Surface reactions on nano- to micrometer sized particles are thought to 
play a key role in the formation of molecules in the ISM Jj. Dust grains in the ISM consist 
of a core of silicates and carbonaceous components and, in dense clouds, these dust particles 
can be covered by icy mantles. The icy mantles contain mainly II2O, but also traces of other 
molecules (e.g., CO, CO2, NH3,CH4, among others) y]. Observed infrared (IR) spectra 
reveal that H2O and CO are the most abundant molecules in the icy mantles in the ISM 



-Iiq|. Recent 



in the ISM 
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ground and space based observations have also detected heavy water (D2O) 
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Ultraviolet (UV) irradiation of an icy grain can photodissociate water molecules and cause 



desorption of the ice. The flux of UV photons in the ISM is low 
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ij] compared with the 



lamp UV photon flux used in the laboratories (which varies between, e.g., (1.1-5.5) x 10 
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photons cm~^s~^ 15|), with photon fluxes of the order of 10^ photons cm~^s~^, which is 
equivalent to roughly one incident photon per month per grain. The photodissociation dy- 
namics is typically computed over a picosecond time scale, and hence, the photodissociation 
by one incident photon is completed before the next photon arrives at the ice. The energies 



of the incident photons ~6-13 eV 

Photodissociation and photodesorption of water in ice are o: 
tronomical observations of gas-phase water in cold clouds 



IGMISI] cover the first absorption bands of water ice. 

interest to understand as- 



19l-|26l|. and also because the 



photoproducts (H and OH) can proceed to react wit 



to the formation of more complex molecules 
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1 co-adsorbed species, which may lead 



27|. In addition, the process is interest- 



ing from a fundamental chemical physics point of view. Most studies of photodissociation 
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processes of molecules on surfaces have focused on (sub) monolayers of species on mostly 
metallic surfaces, not on the thick (~100 monolayer (ML)) ices found in interstellar space. 
Very different processes can occur in this case. 



Several experiments on UV irradiation of amorphous and crystalline H2O ice 



and D2O ice 
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38 



15 
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39| have been carried out using different analysis techniques and 



different light sources. 

In order to obtain insight into the basic molecular processes, the photodissociation of 
H2O molecules in amorphous and crystalline ice in the temperature range 10-90 K has 
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43| . The most important 



been studied using molecular dynamics (MD) simulations 
photodesorption mechanism after photodissociation of water in the top three MLs of the 
ice surface is H atom photodesorption, followed by OH radical photodesorption, and H2O 
molecule photodesorption. The calculated H2O photodesorption probability is due to two 
mechanisms. (1) The direct mechanism: H and OH recombine after HoO photodissociation 
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42|: an energetic 



to form H2O, which eventually desorbs. (2) The kick-out mechanism 
H atom released after photodissociation kicks out one of the surrounding water molecules 
by a transfer of momentum. 

Many key experiments studying water ice photodesorption have been performed for D2O 
rather than H2O. In contrast, the MD simulations have so far been carried out only for H2O. 
In order to better compare with experiments and to identify isotope effects on the photodes- 
orption processes, we present here results of MD simulations of the UV photodissociation of 
amorphous D2O ice at different temperatures. 

In Sec. HI] we present the methods used in this study, in Sec. IIIII the main results in 
comparison with previous UV photodissociation of amorphous H2O ice results are presented, 
and in Sec. |IV]the concluding remarks are given. 



II. METHODS 



A. Potentials 

The total analytical potential energy surface (PES) for the ice in the photodissociation 
calculations is the same as in our previous studies and can be written as follows: 

Vtot = Vice + ^20*-ice + ^20* (l) 



The first term describes the intermolecular interactions between the H2O molecules inside 
the ice excluding the H2O molecule that is photoexcited. These interactions are described 



by the TIP4P potential 4J] with all molecules kept rigid. 

The second term refers to the intermolecular interactions of the photoexcited molecule, 
which is treated as fully flexible, with the rigid ice molecules and the third term is the 
intramolecular potential of the photoexcited molecule. These potential terms also cover all 
interactions involved in the dissociation and possible recombination of the excited molecule. 
The potentials are exactly the same as previously used for photodissociation of H2O ice. All 
details of the potentials and the functions used to switch between different potentials are 



given in Ref. 
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43|. 



B. Amorphous ice surface 

To study the UV photodissociation of D2O ice, we simply changed the mass of the H 
atom to that of the D atom. (Of course, all the interactions that take place in and on the 
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ice during the photodissociation are described with the same potentials employed for the 
UV photodissociation of H2O ice (Eq. [1])). 



Crystalline and amorphous D2O ice surfaces were constructed using the MD method 45 1, 



40|-|43|. 



and using the same procedure and cell parameters employed before to model H2O ice 

Starting from the normal hexagonal ice (Ih) crystalline ice configuration (containing 8 
bilayers (BLs) (16 MLs) with 60 (30) molecules in each ML), the arnorphous ice surface was 



4l|-|43l|. Since the resulting amorphous ice surface has a 



set up at 10, 20, 60, or 90 K using the 'fast quenching' method 
be found in our previous studies 
more irregular bonding structure t 
to MLs is not straightforward 



46l448l|. Further details can 
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43] assigning molecules 



431] a new definition of ML 



lan the crystalline ice surface__ 
43"]. In our most recent sudy 
(binning method 2) was tested and shown to be a more realistic way to assign molecules to 
MLs. This binning (method 2) is used in this study and it consists of choosing a molecule 
and finding the first 23 closest molecules in terms of (x, y) coordinates. This leads to 24 
molecules that are divided in 12 bins of two molecules each, based on their z center of mass 
coordinates. The 12 bins represent the top 12 MLs of the ice in which the molecules are 
allowed to move. 



C. Initial conditions and dynamics 



For each of the top four MLs, all the molecules were chosen to be photodissociated 
and for each molecule 200 different initial configurations were considered. To initialize 



the trajectories 



50l | fitted to the ground- 



49|, a Wigner phase-space distribution function 
state vibrational wavefunction of gas-phase water is used. In the case of heavy water the 
trajectories are initialized by using the corresponding Wigner distribution of gas-phase heavy 



water, which has the same functional form as for gas-phase water (Eq. 5.13 in Ref. 



50j), but 



with aDjO = aHaO (Eq.5.16 in Ref. |50|]). The initial coordinates and momenta of the 
atoms from the dissociating molecule are sampled using a Monte Carlo procedure. Then, 
a Franck-Condon excitation is performed and the system is put on the first electronically 



excited state, on the DK A^Bi PES 
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53|. 



The excitation energies are calculated by computing the energy difference between a D2O 
ice with an excited molecule and one with a ground state D2O molecule (both molecules 
with the same coordinates). The calculated ID^O amorphous ice spectrum is shifted 0.02 eV 



with respect to that for H2O amorphous ice 



4l| 



To simulate the dynamics of a photodissociation event, Newtons's equations of motion 
are integrated in time with a time step of 0.02 fs and a maximum time of 20 ps. The 
stop criterion and the six final outcomes after UV photodissociation of D2O are analogous 
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43|: (1) desorption of D while 



to those described previously for H2O photodissociation 
OD is trapped inside or on the ice, (2) desorption of OD while D is trapped inside or on 
the ice, (3) desorption of both D and OD, (4) D and OD are both trapped inside or on 
the ice, (5) D and OD recombine and form a D2O molecule which either desorbs or (6) is 
trapped inside or on the ice. Besides these six outcomes, an additional channel is possible 
where D2O desorbs through the so-called 'kick-out' mechanism 
a molecule desorbs from the ice by momentum transfer from an energetic D atom resulting 
from photodissociation of a neighbouring photoexcited molecule. 

We calculate the probabilities of the outcomes per absorbed UV photon in a specific ML 
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42|. This occurs when 



i and its standard errors (e'=yP' ■ (1 — P^)/N, where is the total number of trajectories 
simulated in ML i) at all ice temperatures for the top four MLs (the error bars in the figures 
and tables correspond to 66% confidence intervals). However, not all of the UV photons 



that arrive at the ice are absorbed in these monolayers. Andersson et. al. 



42| estimated the 



absorption probability per ML {P^l') to be about 7x10 ^ using an absorption cross section 



42| ) . In the case 



of about 6x 10~^®cm~^ (for more details we refer to the Appendix A in Ref. 
of heavy water amorphous ice we have assumed the same -P^^. The total photodesorption 
yield (Y) can be calculated from the calculated photodesorption probabilities per absorbed 
UV photon in a specific ML i, by multiplying this probability P^^^ with the probability that 
the photon makes it to ML i and the probability that the photon is absorbed in a given 
ML (-P^^), and summing the resulting yields per ML over the considered MLs. This is 



summarized in the following equation 



43|: 



ML 
abs 



(2) 



1=1 



III. RESULTS AND DISCUSSION 



A. D atom photodesorption 



The probabilities of all the different outcomes following photoexcitation of one molecule 
in the ice have been calculated, but we only report those concerning photodesorption (i.e., 
outcomes 1, 2, 3, and 5) and for the top 4 MLs because, according to previous MD simu- 



lations at different ice temperatures 
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42| . the photodesorption mainly takes place after 
photoexcitation in these monolayers. Thus, the outcome probabilities strongly depend on 
the monolayer in which the photoexcited molecule is initially located: the photoexcitation in 
ihe top MLs leads mainly to photodesorption, while deeper into the ice it leads to trapping 
4o|-|43|. 



The average Datom photodesorptiou probabil.ty and the average H atom photodesorp- 

tion probability [43] taken over the top four monolayers (e.g., <PHdes>=Z]f=i -pHdes/^) are 
plotted in Fig. [TJ The average H photodesorption probability is somewhat larger than that 
of D, by about 6 % for all Tice studied. The total deuterium (hydrogen) atom photodesorp- 
tion probability is calculated by summing over two different processes: one in which the D 
(H) atom desorbs while the OD (OH) stays trapped in the surface, and one in which both 
photofragments desorb from the ice surface (outcomes 1 and 3). 

The deuterium and hydrogen atom photodesorption probabilities are higher than the OD 
(OH) and D2O (H2O) photodesorption probabilities in the uppermost monolayers of the 
D2O (H2O) amorphous ice (see Fig. [2] where the H and D photodesorption probabilities are 
displayed versus ice temperature and ML). This is because D and H atoms are smaller and 



are formed with higher kinetic energies immediate 
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y after D2O and H2O photodissociation. 



-l43|. 



which facilitates the desorption of these atoms 

The average photodesorption probability of D is smaller than that of H mostly because 
the probabilities of D atom photodesorption in the third and fourth monolayers (Fig. [2]) 
are smaller than those for the H atom. This trend is expected because D is heavier than 
H. Therefore, the efficiency of energy transfer between D and D2O molecules is larger than 
the corresponding efficiency between H and H2O molecules. If the photoexcited molecule is 
isolated (i.e., in the absence of the surrounding ice), the initial kinetic energy in which D 
and H atoms are formed after D2O and H2O photodissociation should be similar in order to 
achieve energy conservation. But, in the presence of ice, D atoms lose more kinetic energy 
than H atoms when they interact with the surrounding molecules due to the larger efficiency 
of energy transfer, and they are therefore less able to penetrate the ice when moving through 
the upper layers. 
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The dependence on ice temperature is negligible. The average D atom photodesorption 
probability (Fig. [1]) is almost constant at ~54 %. However, this probability depends on the 
ML where the photoexcited molecule was initially located (Fig. In the top two MLs 
the probability is high (~90 % to ~70 %), but it drops in the third ML and further below, 
because other processes such as trapping are in competition. Trapping becomes important 
because deeper in the ice the structure is more closed and the molecules from the ice above 
can impede the D atom from reaching the ice surface. 

D atoms travel through the D2O ice at 90 K by an average distance of 8.4 A before they 
become trapped, whereas H atoms travel around 9.1 A. The OD and OH radicals travel 
2.2 A and 1.9 A, respectively. The recombined D2O and H2O move on average a distance of 
1.8 A and 2.0 A, respectively. Therefore, the mobility of the photofragments inside the ice 
is slightly affected by the mass of the photofragments: H atoms move further than D atoms 
until they become trapped, because H is lighter than D. The maximum distances travelled, 
which are about tens of angstroms should enable reaction with other species trapped in the 
ice. This could explain the formation of more complex molecules in the ISM. 

B. OD radical photodesorption 

The second main photodesorption mechanism in the uppermost MLs of the ice is OD 
photodesorption. Fig. [3] shows that the average of the OD photodesorption probabilities 
taken over the top 4 MLs is larger than that for OH for all ice temperatures studied. To 
calculate the OD (OH) photodesorption probabilities we have summed over the probabilities 
of two pathways: the probability of the channel in which the OD (OH) radical desorbs while 
the D (H) atom remains trapped in the ice surface, and the probability of the channel in 
which both photofragments leave the ice surface (outcome 2 and 3, respectively). 

11 
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FIG. 1: The probability of D atom (dashed hne) and H atom (sohd hne) photodesorption averaged 
over the top fom' MLs per absorbed UV photon is shown as a function of ice temperature. H atom 
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FIG. 2: Total probability of D atom (dashed line) and H atom (solid line) 
(per absorbed UV photon) versus temperature, for the uppermost four MLs. 



43l | photodesorption 
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The probabilities per monolayer are plotted in Fig. H] versus ice temperature. Those for 



OD are much larger than those for OH 43|] in the top two monolayers, which gives rise 



to a larger total average OD photodesorption probability (Fig. E]). In the absence of the 
surrounding ice, the water fragments (D and OD, H and OH) have to obey momentum 
conservation {px=—Pox, -^=H, D) and energy conservation {Ex+Eox=Ex20=E , E being 
the initial available energy ii^exc— -Ediss(X20), X=H, D; the excitation energy Eexc is in the 
range 7.5-9.5 eV with a peak at 8.6 eV, and the dissociation energy £'diss(X20) ~ 5.4 eV 



53|), leading to the following equation: 



^ 2 ^ /Q\ 



Thus, if the molecule is isolated and dissociates, the OD radicals will be formed with 
higher initial translational energy than the OH radicals (Eod=-E'/10 and Eoh=-£'/18, which 
leads to £"00 ~ I-Sx-E'oh) according to Eq. |3l Because the OD radicals have a higher initial 
translational energy than the OH radicals in the uppermost monolayers, they leave the ice 
surface more easily. 

The OD photodesorption probability decreases with increasing depth, similar to the OH 



photodesorption probabihty 
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43| . In the third and fourth MLs, the OD and OH pho- 
todesorption probabilities drop to less than Ix 10~^ because the OD and OH radicals do not 
have enough translational energy to escape from the ice surface. 

An oscillatory effect is observed when the OD and OH photodesorption probabilities are 
plotted in MLs 1-4 versus ice temperature (Fig. H]). We attributed these oscillations in our 



previous paper 43|] to the irregular nature of the amorphous ice surface, which makes it very 
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complicated to assign molecules to specific MLs, and also to the finite sample size of about 
30 molecules per ML. 

The average OD and OH photodesorption probabilities increase with ice temperature by 
~24 % and ~25 %, respectively, from 10 to 90 K. If longer time scales in our simulations 
could be considered, a stronger dependence on ice temperature would be expected, because 



processes like thermal diffusion and thermal desorption are more efficient at higher 
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FIG. 3: The probability of OD radical (dashed line) and OH radical (solid line) photodesorption 
averaged over the top four MLs per absorbed UV photon is shown as a function of ice temperature. 



OH radical results from Ref. 



43] 



C. D2O molecule photodesorption 



1. Kick-out vs. direct mechanism 



The third photodesorption channel upon UV photodissociation of D2O amorphous ice is 
D2O molecule photodesorption. This is due to two mechanisms: the direct and the kick-out 
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FIG. 4: Total probability of OD radical (dashed line) and OH radical (solid line) 43] photodes- 
orption (per absorbed UV photon) versus ice temperature, for the uppermost four MLs. 

mechanism. The direct mechanism consists of the recombination of the D atom and OD 
radical that leads to the formation of an energetic D2O molecule that eventually desorbs 
(outcome 5). In this situation the D2O molecule has a high probability of desorbing from 



the ice in a vibrationally excited state 42|. The kick-out mechanism takes place after the 
photodissociation of a D2O molecule when an energetic D atom transfers momentum to 
one of the surrounding D2O molecules, which is then likely to desorb vibrationally cold 
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43| . Since in our model we cannot quantify the energy transfer from X to internal 
modes of X2O, because the kicked out molecule is treated as internally rigid, we have carried 
out quasi-classical trajectory calculations on the isolated X-X2O system at different incident 
X atom kinetic energies (X=H, D). Here we only report the results for Ex=l-5 eV, because 
it is the average kinetic energy with which H and D atoms kick a surrounding molecule out of 
the ice at TirB=10 K. Further details on these calculations will be reported in the paper we are 
preparing 



54j |. Gas phase collisions at Ex=l-5 eV lead to final average vibrational energies 
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of ~0.16 eV for H2O, and ~0.77 eV for D2O, whereas in the ice system recombined H2O 
and D2O molecules desorb from the ice surface with average ro-vibrational energies of 5.3 
and 5.4 eV, respectively. Thus, the kicked out molecules are much more likely to be formed 
in states with lower vibrational energies than the molecules formed after recombination of 



the photofragments. Yabushita et al. 



observed photodesorbed 



29[, and Hama et al. 

X2O in the ground vibrational state by using resonance-enhanced multiphoton ionization 
(REMPI) detection methods and classified them as kicked out. They speculated that when 
the X atom kicks out a X2O molecule, most of the energy is transferred into translation and 
much less into internal energy because the X atom hits the X2O molecule close to the center 
of mass. This is supported by our simulations (e.g., ~89 % of the trajectories classified as 
kicked out occurs when the X atom kicks the X2O molecule close to the oxygen atom at 
Tice=10 K) (X=H, D). Isotope effects are also observed in our calculations: the efficiency of 
energy transfer from D to intramolecular vibrational modes of D2O molecule in gas phase 
collisions is larger than the corresponding efficiency for collisions of H with H2O, because 
D is heavier than H and because the vibrational frequencies of D2O are lower than those of 
H2O. 

Fig. El^a) shows the average of the D2O photodesorption probabilities compared with 



those for H2O 



43| over the top four monolayers versus ice temperature. The average values 



for D2O and H2O [43|] due to the kick-out mechanism are displayed in Fig. |5]^b), and those 
due to the direct mechanism in Fig. El^c). 

The total average D2O photodesorption probabilities are larger than those for H2O 
(Fig. |5]^a)) because of the large contribution of the D2O kick-out mechanism (Fig. ISt^b)). 
For every ice temperature studied, the D2O photodesorption probability due to the kick-out 
mechanism is larger than that for H2O (Fig. |5]^b)), because in the case of heavy water, the 
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D atom formed after D2O photodissociation has a larger momentum (by about a factor 
a/2) than the H atom formed after H2O photodissociation. In addition, the efficiency of 
momentum transfer from D to D2O is larger than that from H to H2O, so that the kick-out 
mechanism is much more successful for D2O. The effect occurs mainly after photoexcitation 
in the second and third monolayers (see Figs. Mjo) and[6](c)), because a D atom produced 



llM. The kick-out 



in these MLs is more likely to kick-out a molecule located above it 
photo desorption probabilities are much lower in the first and fourth MLs (Figs. El^a) and 
ini^d)), where there are no important differences between D2O and H2O. 

We have calculated the probabilities of the parallel outcomes that take place in coinci- 
dence with the kicking out of a D2O molecule by an energetic D atom. These probabilities 
are summarized in Table [H We do not observe differences with those calculated for the H2O 
kicked out molecules from H2O ice (Table I in Ref. 43|). The most dominant simultaneous 



process is that where the D atom that kicks out the D2O molecule also desorbs, while the 
OD fragment is trapped (Table [T]). The next most important parallel processes are those in 
which the photofragments recombine and form a D2O molecule that remains trapped in the 
ice, and those in which both photofragments are trapped inside the ice at separate locations. 
It is also possible that two molecules desorb at the same time, i.e., the kicked out molecule 
and a recombined molecule, but this process occurs with very low probability. 

2. Trends with ice temperature 

The total D2O photodesorption probability increases faster with ice temperature than 



43| . going from 10 to 90 K. The average 



that for H2O (Fig. EJ^a)): by 130 % vs 30 % 
desorption probability due to the direct mechanism (Fig. Et^c)) is relatively small, and there 
are no differences between the average D2O and H2O photodesorption probabilities: both 
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TABLE I: Probabilities averaged over the top four monolayers of the outcomes that take place in 
coincidence with the kicking out of a D2O molecule for each ice temperature. Overall probabilities 
can be obtained by multiplying the probabilities shown with the probabilities for the kick-out 
mechanism, see Fig. [5{b). 



Tice / K 
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0.257 ± 0.041 


0.257 ± 0.041 
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0.412 ± 0.040 


6.8 ± 6.7 


6.8 ± 6.7 


0.108 ± 0.026 


0.459 ± 0.041 


6.8 ± 6.7 


60 


0.511 ± 0.031 








0.211 ± 0.025 


0.278 ± 0.027 





90 


0.587 ± 0.028 


3.2 ± 3.2 


3.2 ± 3.2 


0.167 ± 0.021 


0.237 ± 0.024 


3.2 ± 3.2 



rise with ice temperature by only ~30 % going from 10 to 90 K. However, the D2O average 
kick-out photodesorption probability increases strongly with ice temperature (by ~180 % 
from 10 to 90 K), whereas the H2O kick-out average photodesorption probability shows 
a much weaker increase with ice temperature (by ~50 % from 10 to 90 K). Thus, the 
stronger trend with ice temperature for D2O results from the increase of the D2O kick-out 
photodesorption probability (Fig. El^b)). This probability increases with ice temperature 
because the molecules have higher initial kinetic energies at higher ice temperatures, which 
promotes the desorption of the surrounding molecules through the kick-out mechanism. 

At higher ice temperatures the D2O kick-out mechanism has a high probability (per ab- 
sorbed UV photon): 2.8 % at 60 K after photoexcitation in the 3rd ML (Fig. EJ^c)) and 2.2 % 
at 90 K after photoexcitation in the 2nd ML (Fig. |6t^b)). The probabilities for D2O and H2O 
desorption through the kick-out mechanism in the top 4 MLs (Fig. |6]) show an oscillatory 



18 



dependence on ice temperature, as also seen before for OD and OH photodesorption proba- 
bilities versus ice temperature and monolayer (Fig. H]). These oscillations are attributed to 



the corrugation of the amorphous ice surface (for more details see Ref. 



43|). 
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FIG. 5: (a) The total D2O and H2O [43] photodesorption probabilities, (b) the D2O and H2O [43| 



photodesorption probabilities due to the kick-out mechanism, and (c) the probabilities due to the 
direct mechanism (per absorbed UV photon) versus ice temperature, all averaged over the top four 



MLs in which the photoexcited molecule resides. 



D. Energies of the kicked out molecules 

The average translational and rotational energies taken over the top 4 MLs of the D2O 



and H2O 



43| kicked out molecules are plotted versus ice temperature in Fig. [71 together 



with the corresponding experimental values for H2O photodesorbed molecules (f =0) at 90 K 



291 . The translational energies tend to increase with ice temperature, because the energy of 
the ice also rises. The final rotational energies are low and do not display any dependence 
on ice temperature. The calculated translational and rotational energies do not show a 
significant isotope effect. Our calculations cannot say anything about the vibrational state 
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FIG. 6: Probabilities of D2O molecule (dashed line) and H2O molecule (solid line) [43[ photodes- 



orption due to the kick-out mechanism upon photoexcitation in (a) the first ML, (b) the second 
ML, (c) the third ML, and (d) the fourth ML, (per absorbed UV photon) versus ice temperature. 

of the kicked out molecules, because the molecules that are not photoexcited are kept rigid 
in our model. However, it seems unlikely that the kicked out molecules would emerge highly 
vibrationally excited. The average translational and rotational energies at 90 K of the H2O 



kicked out molecules are 0.29 and 0.044 eV, respectively [43[, and of the D2O kicked out 
molecules 0.27 and 0.021 eV, respectively. These results are in good agreement with the 
experimental translational and rotational energies of the H^O desorbed molecules in their 



ground vibrational state as measured by Yabushita et al. 29|] at 90 K (0.31 and 0.039 eV, 
respectively), and also with the experimental translational and rotational energies measured 
by Hama et al. [38|] for H2O and D2O ices at 90 K (0.31 and 0.047 eV, respectively). Hama 



et al. 



38 



did not observe differences between the energies of desorbed H2O (f=0) and 



desorbed D2O (f =0) at 90 K, in agreement with our calculations. 
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FIG. 7: Calculated average translational and rotational energies of the kicked out D2O molecules 



versus ice temperature, of the kicked out H2O molecules 43|, and experimental average trans 



lational and rotational energies of H2O molecules desorbed in their ground vibrational state at 



E. Total (OD + D2O) photodesorption yield and comparison w^ith experiments 



The calculated average of the total (OX + X2O, for X=H or D) photodesorption probabil- 
ity per absorbed UV photon is larger for D2O ice than for H2O ice at every ice temperature 
(Fig.i). 



In our previous study 



43'], we compared the total photodesorption yield of the O con- 



taining species (H2C)des + OHdes) with the total experimental photodesorption yield (Eq. 4 
in Oberg et al. 15]). In the calculation of this yield, the deuterium (or hydrogen) atom 
photodesorption is not included (outcomes 2, 3, 5, and the kick-out), because D (or H) 
was not detected in the experiments by Oberg et al. jl5|. The photodesorption yields (i.e., 
photodesorption probabilities per incident photon) have been calculated through Eq. [2] (see 
Sec.Hm. 
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FIG. 8: The probability of (OD^es + D20des) photodesorption (dashed hne) and (OHjes + H20des) 



photodesorption (sohd hne) averaged over the top four MLs per absorbed UV photon is shown as 
a function of ice temperature. The (OHdes + H20dcs) photodesorption results are taken from 



Ref. 



43| 



Although the actual experimental photodesorption yields (Table [TTl) 



ice, Oberg et al. 



15| were for D2O 



15| applied the results for D2O ice to H2O ice (Eq. 4 in Ref. 



15|). because 



they found that the total photodesorption yields from D2O and H2O were indistinguishable 
(i.e., no isotope effects) within the experimental uncertainties (60%') at 18 and 100 K. 



Table |TT] contains the total experimental photodesorption yield 



(H20dcs + OHdes) 



15| . the computed total 



43| . and the computed total (D20des + ODdes) photodesorption yield 



per incident photon, the ratio ^ between the experimental yield and the calculated (H20des 
+ OHdos) yield, and the same for the (D20dos + ODdes) yield, at all ice temperatures. ^ 
increases from 3.0 to 5.9 for H2O ice and from 1.3 to 2.3 for D2O ice. Thus, the computed 
(D20dcs + ODdes) photodesorption yield is significantly larger than the computed (H20des 
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+ OHdcs) yield at all ice temperatures, as also illustrated in Fig. [HI and compares better 
with the experimental D2O photodesorption yield than the yield calculated for amorphous 
H2O ice [43]. 

The agreement between theory and experiment is better at low ice temperatures, (e.g., 
,^=1.3 at 10 K and ^=2.3 at 90 K). This trend strengthens our explanation that the dif- 
ference between theory and experiments can be due to long time scale processes promoted 
by prolonged irradiation effects leading to an accumulation of radicals, thermal desorption 
and thermal diffusion. At higher ice temperatures, the photofragments, possibly formed by 
different photodissociation events, become more mobile, allowing them to recombine and 
eventually desorb as a consequence of the excess energy. Some of the OD photofragments 
that are trapped deeper in the ice, could probably desorb at long time scales due to a higher 
diffusion rate at higher Ti^e- In our simulations we can only reach the picosecond time scale. 



therefore these kind of secondary processes are beyond the scope of our simulations 43 1 
Another difference is the UV wavelength covered by the lamp used in the experiments 15 1 
This UV lamp includes Lyman-a photons which can excite H2O to the B state whereas our 
calculations consider only the A state 



4l| . Given the experimental uncertainties and our 
approximations (such as the use of a gas phase PESs for the H2O intramolecular interactions, 
the freezing of the intramolecular degrees of freedom of the surrounding molecules, and the 



short time scale of our simulations 



41 



43l|). the experimental and calculated probabilities 
may be considered to be in reasonable agreement. An important result for astrochemists 
is that the computational results fall within the range of the photodesorption probability 



per incident photon (1x10 ^-3.5x10" 



55N59| used to model astrophysical environments. 



However, our calculations suggest that the computed total photodesorption yield can be 
different for H2O and D2O ice, in contrast to the experimental results. 
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TABLE II: Experimental ^, theoretical (OHdes + H20des) |43|, and theoretical (ODdes + D20des) 
photodesorption yields per incident photon, the experimental yield/theoretical yield ^ (H2O) and 
the experimental yield/theoretical yield ^ (D2O) at all ice temperatures. 



Tice / K 



Exp. (OXdes + XaOdes), X=H or D 
xlO"^ 



(OHdes + H20de 



(ODdes + DaOde 



e (H20) 



c (D20) 



10 
20 
30 
60 
90 



1.62 ± 0.48 
1.94 ± 0.56 
2.26 ± 0.64 
3.22 ± 0.88 
4.18 ± 1.1 



0.54 ± 0.04 
0.57 ± 0.04 
0.71 ± 0.09 



0.71 ± 0.05 



1.27 ± 0.063 



1.46 ± 0.066 



1.57 ± 0.068 



1.83 ± 0.074 



3.0 



3.4 



3.2 



5.9 



1.3 



1.3 



2.1 



2.3 



Calculated from the empirical fit of the total photodesorption yield, Eq. 4 in Oberg et al. 



15j. 



IV. CONCLUSIONS 

In this work we have studied the processes following UV photodissociation of D2O in 
amorphous heavy water ice and compared them with previous UV photodissociation results 
in H2O ice at different ice temperatures, in order to investigate isotope effects in photodes- 
orption. 

D atom photodesorption is the most important desorption process in the uppermost MLs 
of the ice, like H atom photodesorption. The average D atom photodesorption probability is 
smaller than that of the H atom, because in the top two MLs of the ice both the H and D atom 
can easily escape from the ice surface, but if the atoms are located in the third and fourth 
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MLs, the D atom is less likely to penetrate the upper ice layers due to more efficient coUision 
energy transfer to D2O. Therefore, the D atom photodesorption probabilities in these MLs 
decrease, and the same is then for the average D atom photodesorption probabihty. The D 
and H atom photodesorption probabilities do not show any dependence on ice temperature. 

OD and OH radical photodesorption constitute the second most important desorption 
channel in the top two MLs of the ice. Deeper into the ice the probabilities of these processes 
drop because OD and OH radicals do not have enough translational energy to desorb from the 
surface. The average OD photodesorption probability is higher than that of OH. This trend 
can be explained by the initial translational energy of OD being higher by about a factor 1.8, 
a result obtained if the photoexcited molecule is considered to be isolated (i.e., in the absence 
of the surrounding ice) and the laws of momentum and energy conservation are applied. The 
average OD photodesorption probabilities increase smoothly with ice temperature, by about 
24 % from 10 to 90 K. 

The third most important desorption mechanism is D2O and H2O photodesorption. This 
process takes place either by direct desorption of the photoexcited molecule after the re- 
combination of D (H) and OD (OH) or by indirect desorption due to an energetic D (H) 
atom which transfers part of its momentum to a surrounding molecule that is kicked out 
from the ice surface. The average photodesorption probability is higher for D2O than for 
H2O at all ice temperatures considered. This trend is due to the contribution of the kick-out 
mechanism, which is much more important for D2O than for H2O. This result is expected 
because the photoproduced D atoms have higher average momentum (by about a factor v^) 
than the H atoms, and because energy transfer in D-D2O collisions is more efficient than 
energy transfer in H-H2O collisions. The kick-out mechanism mainly takes place when the 
photoexcited molecule is initially located in the second and third MLs of the ice. Photodis- 
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sociation leads to an energetic D (H) atom that can transfer its momentum to a molecule 
located above it, which will desorb from the ice if it has enough kinetic energy. The average 
direct photodesorption probability (which involves recombination) does not show any iso- 
tope effect. The average total D2O photodesorption probability tends to increase with ice 
temperature faster than that of H2O: by ~130 % vs ~30 %, from 10 to 90 K. 

Experiments show, and a consideration of the mechanism suggests, that the kicked out 
molecules leave the surface vibrationally cold. In contrast, the molecules that desorb due to 
the direct mechanism are formed vibrationally excited. The average translational and rota- 
tional energies in which the D2O and H2O molecules desorb due to the kick-out mechanism 
have been calculated and compared with the corresponding experimental values at 90 K. 
The agreement between our MD calculations and the experimental measurements is good, 
and leads to the conclusion that the final energies with which the kicked out molecules are 
formed do not display an isotope effect. 

We have also estimated the total photodesorption probabihty (ODdes + D20des) per 
incident photon from the total photodesorption probabilities per absorbed UV photon, and 
compared this quantity with the previously calculated values for (OHdes + H20(ies), and with 
the available experimental yields. Our total photodesorption probability for D2O compares 
better with the experimental photodesorption yield than that for H2O, and also better at low 
ice temperatures. Presumably at higher ice temperatures long time scale processes become 
increasingly important, such as diffusion and thermal desorption, which are not covered in 
our picosecond simulations. 

Current experiments cannot distinguish between (OH + H2O) and (OD + D2O) yields 
within the experimental uncertainties of 60%. More accurate future experiments may reveal 
the isotope effects predicted here. 
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